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Modelling the composition of a young star cluster ejecta 
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ABSTRACT 

We have computed with a fine time grid the evolution of the elemental abundances of He, 
C, N and O ejected by young (f < 20Myr) and massive (M = 10^ M©) coeval stellar clus- 
ter with a Salpeter initial mass function (IMF) over a wide range of initial abundances. Our 
computations incorporate the mass loss from massive stars (M > 30 Mq) during their wind 
phase including the Wolf-Rayet phase and the ejecta from the core collapse supernovae. We 
find that during the Wolf-Rayet phase (f < 5 Myr) the cluster ejecta composition suddenly 
becomes vastly over-abundant in N for all initial abundances and in He, C, and O for initial 
abundances higher than l/5th Solar The C and O abundance in the cluster ejecta can reach 
over 50 times the solar value with important consequences for the chemical and hydrodynam- 
ical evolution of the surrounding ISM. 

Key words: galaxies: abundances - stars: abundances - stars: mass-loss - stars: Wolf-Rayet 
-stars: supernova - galaxies: stars clusters 



1 INTRODUCTION 

Massive young stellar clusters are ideal laboratories for research 
into the evolution of massive stars and their interaction with 
their surrounding interstellar medium (ISM). These luminous and 
rapidly evolving massive stars supply most of the young cluster 
UV radiation that creates the encompassing Hii region and a large 
amount of mass and mechanical energy in the forms of supernova 
ejecta and stellar winds particularly during the WR phase. These 
massive starforming regions can eject durin g their first 10 Myr of 
evolution about 20% of their initial mass jLeithereretaill 19991) 
(hereinafter STB99) mostly in the form of newly synthesized C, N 
and O that by mass represent most of the heavy elements. 

Massive young stellar clusters are ubiquitous particularly 
among late type galaxies. Their stellar wind phase can result in 
a supergalactic wind affecting the nearby intergalactic medium 
(IGM). It is open to question how the ejected freshly synthe- 
sized heavy elements cools and mix with the ISM, and how long 
this process governed by the cooling time scale, may last. The 
time scale for cooling is strongly dependent on the gas cool- 
ing rate that in turn is dependent on the gas chemical compo- 
sition and density (e.g . Tenorio-Tagle, 1996; Kobulnicky et al., 
1 9971: 'Kobulnicky & Skillman J 19981 : Ivan Ze e & Havneg, l200d : 
Lopez-Sanchez & Esteban. 20101) . Thus evaluating the composi- 
tion of the cluster ejecta and its time evolution is a necessary prior 
step for estimates of the evolution of the cooling function and 
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the computation of cooling and feedback time scales. Although 
work like that of Silich et al. (2001.) underlined the large influence 
that the enrichment of the stellar ejecta can have on the radiative 
cooling of starburst su perbubbles, many researchers are still using 
[Raymond etal] ( Il976t) radiative cooling coefficient calculation for 
Solar abundances when modelling the interaction between the stel- 
lar ejecta and their surrounding medium. 

In this paper we present a set of models designed to calculate 
in detail the first 20 Myr of the evolution and chemical composition 
of a star cluster ejecta on very short time scales, i.e. much shorter 
than the Hii region lifetime and in particular to resolve the WR wind 
phase 

Galactic chemical evolution models traditionally assume that 
the elements ejected by the stellar cluster are incorporated to the 
ISM whe n the corresponding stars die, i.e. at a time equal to their 
lifetimes jPortinari. Chiosi, & Bressanl 1 19981 : iDray & Toull l2003l : 
iDrav et all |2003|) . The shortest time step is usually defined by 
the mean-lifetime of the most massive star. Since this is typically 
around 100-120 Mq, chemical evolution calculations begin nor- 
mally at around 3-5 Myr with comparable time-steps. This way im- 
portant phases of the wind evolution, occurring before 3 Myr or 
short lived like the WR phase are lost or diluted. 

To compute the composition of the ejecta, most chemical 
evolution models use the total yields of el ements due to su- 
pemova explosions, such as those given by Wooslev & Weave^ 
(1995,W W95) or other more recen t works ( Umeda & Nomoto , 
20021: IRa uscher et al.'. 2002; Limongi & Chieffi, 2003; 



Chieffi & Lim ongi. 2004; Frohlich et al., 2006; Heger &_ Woosley 
20 lot) to calculate the total change of the elemental abun- 
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dances. Other computations include the elements 
during the wind phase of massive stars (Maede i 
Mevnet & Maeder, 200 j. iHirschi. Mevnet. & Maedei 



Hirschii i2007i; iKobavashi et al 



200d) . 



ejected 



1992; 



2005 



Table 1. Chai'acteristics of the stellar models 



only the models of 

Portinari. Chiosi, & Bressanl(ll998l hereinafter PCB98) include the 



evolution of both phases, i.e. the yields of co re collapse supernova 
explosions from Woosley & Weaveil ( Il995h and the stellar wind 
yields produced during the evolution of each star. However, PCB98 
computations, as most chemical evolution models, were performed 
with time steps much longer than the lifetime of an Hn region 
therefore missing short lived stages like the WR phase. Moreover, 
the evolution of a supernova progenitor that loses part or most of its 
mass is not the same as a normal main sequence massive star Since 
mass and structure are substantially different in he time of the 
supernova explosion, the associated yields w ill also differ, such as 
it is ex plained in I Woosley et all ( ll993l) : IWooslev. Langer & Weavej 
( I1995L hereinafter WLW93 and WLW95, respectively). To take 
this into account, PCB98 linked the final stage of the star after 
lost mass with the supernova yields though the CO core mass. 
However, they use the WW95 yields instead WLW93/WLW95 
yields. 

The organization of this paper is as follows: In Section 2 we 
use O and WR winds computations to estimate the evolution of the 
ejecta of a young massive stellar cluster In section 3 we compute 
the contribution of the explosive nucleosynthesis. Section 4 gives 
the complete evolution of elemental abundances within the cluster 
and discusses the impact of each contribution phase over the final 
ejecta. Our conclusions are given in Section 5. 



2 STELLAR WIND COMPOSITION 

To compute the evolution of the rate of ejection of elements 
due to stellar winds of a Single Stellar Population (SSP) of 
a given initial metallicity and initial mass function (IMF), we 
have used the tables result i ng from the isochrones calcu lation by 
iBressan. BerteUi & Chio"sil ( Il993h ; iFagotto et aO \\99Aa.h). These 
authors give the amount of mass lost by massive stars during their 
evolution for 7 initial stellar masses: 12, 15, 20, 30, 40, 60 and 
IOOMq. For each initial mass and time step, the tables provide the 
present mass (in solar mass units), the rate of mass loss in Mq yr"' 
(logarithmic scale) and the abundances of the stellar surface for 
H, *He, '"C, '■*N, and '^O. These tables are provided for 6 initial 
metallicities: Z=0,0001, 0.0004, 0.004, 0.008, 0.02 and 0.05. 

In Fig. [T] we plot the evolution of the mass loss rate for all 
stellar masses. The mass loss rate depends strongly on the initial 
stellar mass and composition. It is clear from the figure that the 
lowest mass stars maintain for a long time a low mass loss rate 
while the most massive ones evolve rapidly ejecting a large part 
of their mass in discrete events. An important consequence is that 
a star of solar metallicity with an initial mass of IOOMq ends its 
life with around 7 Mq, while the evolution of star of 12 Mq may be 
followed for almost 20 Myr at a very low mass loss rate that implies 
that its total mass remains roughly constant. 

The mass loss rate is also dependent on the metallicity through 
the semi-empirical relation incl uded in the stellar models (see 
iBressan. BerteUi & Chiosill 19931 for details). We see these differ- 
ences in Fig.[T]where the evolution for the 6 given metallicities stars 
are shown. The lower the metallicity, the smaller the mass loss rate 
and smoother the behavior shown on the mass loss rate evolution. 

In Table [T] we summarize some characteristics of the stellar 
input models: For each metallicity Z, column 1, and initial stellar 
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mass m,, column 2, the final mass m^nd is listed in column 3, and the 
final evolutionary time te„d in Myr is in column 4. Column 5 shows 
the final surface abundance (in mass fraction) of H, Xf ^ (which we 
will use to determine the mass of the Helium core in section 3). 

Since each stellar mass table has a different time range and 
step, we have performed a linear interpolation to obtain values at 
the same times for all masses. The noiTnalized table is provided in 
electronic format for all calculated metallicities and masses. Table|2] 
shows as an example, a few time steps of the most massive star of 
Z=0.02. It gives for the metallicity of column 1, and for the seven 
stellar masses, defined by column 2, the initial mass m, in solar 
mass units in column 3, the time in Myr units in column 4, the 
mass loss rate, m, in units of solar mass per year, in column 5, the 
elemental abundances of H, "He, '^C, '''N, and "^oQ as fractions in 



' For the sake of simplicity we write He, C, N and O for ''He, '^C, 
and "O along the text. 
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Figure 1. Evolution of the mass loss rate for the stai's given by the Padova group for different metallicities Z as labelled. Color and type of lines mean a 
different value of stellar mass as labelled. 



Table 2. Surface abundances of the stellar models for each mass for Z=0.02. Similar complete tables for the six metallicities and seven stellar masses will be 
provided in electronic format 
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Figure 2. Input stellar parameters for Z=0.02. The evolution of the stellar mass is shown in panel (a) and the evolution of the surface abundances is shown 
in panels (b) to (f) for the same initial stellar masses as labelled in Fig. 1 . The dots represent the actual results from the stellar evolutionary models while the 
lines are the interpolation in time used in this work. 



mass, in columns 6 to 10, and the stellar mass m(t) at a given time 
t, in column 1 1 . 

The evolution of the stellar masses and surface abundances 
is shown with dots in Fig.|2] Each type and color of line indicate a 
stellar mass as labelled. The lines shows the results of the numerical 
interpolation used in the following sections. In panel a) we see how 
drastically the stellar mass decreases when m > 30 Mq in times 
as short as 5 Myr. The evolution of the stellar surface abundances 
for H, He in total mass fraction and C, N, and O, as abundances in 
mass, X, is shown in panels b) to f) of the same Fig. |2] The sur- 
face abundances of C, N and O show a large increase following the 



start of the stellar winds revealing the the product of first hydrogen 
and then helium burning. This fact combined with a depletion of H 
in the ejecta means that, if these abundances were represented as 
abundances in number, 12 + log (X/H), they would be very high. 

To compute the evolution of a stellar cluster ejecta we have 
assumed that the cluster stellar mix consists of a coeval popula- 
tion or single stellar population (SSP) where all stars were created 
simultaneously and with the same metallicity. By using the normal- 
ized tables of the previous paragraph, it is easy to calculate in each 
time step the contribution of each star, m, weighted by the number 
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Figure 3. Evolution of the instantaneous abundances in mass for a stellar cluster due to stellar wind integrated for a Salpeter IMF: a) He, b) C, c) N and d) O for 
Z=0,05, 0.02, 0.008, 0.004,0.0004 and 0.0001 with different coded lines as labelled in panel a. In each panel the dotted grey line marks the initial abundance 
for Z=0.02 and the arrows the initial values for the other metallicities 



of stars in its mass range, given by the initial mass function (I)(m). 
Thus, for each element i and each time t: 



where 



Jnmup 
mlow ^ Ar 



(m, t')^(m)dmdt' 



e,_,(m, /') = XSi(m, f)m(f) 



(1) 



(2) 



XSi(m, t) is the surface abundance of each element ; and m{t) 
is the mass loss rate for each stellar mass m in every time t. 

We have performed the calculations for 6 different initial 
mass functions: 1) A Saloeter (1955) law, <I)(m) oc m"^, with an 
exponent x = -2.35 (hereinaft er SA L), and 5 others from: 2) 
iMiller & Scalol (hereinafte r MIL|1979|) . 3) iFerrini. Penco & Pallal 
(ll990Lher einafter PER), 4) iKroupd (hereinafter KRO l200lh . andS) 
IChabriei] (hereinafter CHA 12003'); all of them with limits m/o„ = 
0.15 and m„,„v = 100 M© and 6) a Salpeter law with limits m/ow = 
1.00 and m„,a> = 100 Mq as the one used in STB99. Our tables are 
therefore calculated for all these IMFs, but in the next figures only 
Salpeter results are shown since it is a widely used IMF. 

The integration is done for the whole mass range of the IMF 
in each time step.To integrate in time we have chosen a time step 
6t = 0.01 Myr (small enough to follow the rapid evolution of the 



mass loss process without losing any phase). To integrate in mass 
it is necessary to be available a grid with a wide range of masses. 
So, we have performed a careful interpolation in mass using the 
7 existing tables. The method to obtain the mass loss rate for any 
mass value it is not straightforward, since it shows abrupt changes 
in small time scales. So we have used the tables of the mass loss 
rate to calculate the actual mass m(t) in any time with high accu- 
racy. Then we have interpolated between these values for obtaining 
a new curve for each mass m, and finally we compute the mass loss 
rate from this m{t). To calculate the elemental abundances, we have 
taken into account the different phases of each element abundance, 
interpolating between two known masses to obtain the points limit- 
ing these phases for each mass m. All computations were done for 
a SSP cluster with a total mass in stars of 10^ Mjf|. 

In Fig. [3] we show the resulting evolution of the instantaneous 
ejecta abundances X for He, C, N and O in mass as before. Initial 



^ It is necessary to take into account that we will give averaged val- 
ues for this stellar cluster mass. For stellar clusters less massive than 
Iff^Me, the stochastic effects over the initial mass function are important 
liCerviiio & Luridianal I2O0I |2006|) and they may change the stellar mass 
distribution and the corresponding results compared with those we obtain. 
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Figure 4. Evolution of abundances in mass of tlie accumulated masses ejected by a stellar cluster for a) He, b) C,c) N and d) O. The metallicities are shown 
with the same color codes as in the previous figure. The dotted line and the arrows have the same meaning than in Fig. 3 



Table 3. Accumulated Masses ejected by a stellar cluster of 10^ M© during its wind phase for 20Myr for Z=0.02 and a Salpeter IMF. The complete tables for 
all times, metallicities and IMFs will be provided in electronic format 
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abundances are the ones included in the original files for Z=0.0001, 
0.0004, 0.00 4, 0.008, 0.02 and .05. If we take the solar abun- 
dances from lAsplund et al. I ( l2009l) as reference, which implies a 
total Zo = 0.0142, the Z=0.02 abundances, indicated in the figure 
as dotted grey lines, are in fact ~ 0. 15 dex higher than the adopted 
Solar values. If we analyze the results for Z=0.02 (red lines) we 
see that cluster ejecta abundances are not always equal to the initial 
value. It is evident that ejecta abundances show strong variations in 



the WR stars phase, that may reach up to two orders of magnitude 
over the initial value for each metallicity. Thus, He and N increase 
in a first phase, when C and O decrease. Then N and He decrease 
again just when C and O increase. 

This efi"ect is strongly dependent on the original metallicity of 
the stellar cluster as we may see by comparing with the other metal- 
licity lines. The initial values for the non-solar cases are indicated 
with arrows around the values, +0.3, -0.4,-0.7, -1.7 and -2.3 dex. 
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Figure 6. Evolution of the instantaneous mass in the higher mass stellar 
models. The horizontal lines give the WLW93/WLW95 final masses. 

For the two models with the lowest metallicity, the ejecta abun- 
dances, in particular C and O, do not differ very much from the 
initial values except for the decrease around 4-5 Myr. However, for 
higher initial metallicities, all lines show a strong sudden increase 
at around 4.5-5 Myr reaching values 1 or 2 orders of magnitude 
larger than the initial value and lasting few million years. 

For N the behavior is however different. For all models the 
N abundance in the ejecta increases at 3 Myr, even at the lowest 
metallicities. After about 5 to 8 Myr of evolution, when the clus- 
ter tumoff mass is below 25 M© the ejecta metallicity asymptoti- 
cally approaches the initial value for all initial compositions, ex- 
cept again for N which maintains a higher value than the initial 
one. Probably this may be explained by the fact that the convective 
envelope dredges up the modified composition of CNO from the 
inner parts of the star up to the surface during the Red Supergiant 
(RSG) RSG phase, which enhances the N mass fraction, sligthly 
reducing C and O abundances. 

The resulting accumulated ejected masses for the cluster in 
every time step are given in Table |3] For each IMF, column 1, and 
metallicity, given in column 2, we give the time step in Myr in 



Table 4. Relation stellar mass-He core mass for Padova stellar tracks 





0.008 


Z 

0.02 


0.05 


m. 




Mh, 




30 




12.05 


9.84 


40 


17.04 


14.50 


15.56 


60 


25.44 


27.80 


20.13 


100 


45.14 


42.20 


28.44 



column 3, the total ejected mass in column 4, and the accumulated 
ejected mass of the different elements in columns 5 to 9 for H, He, 
C, N, and O. We show an example here for solar metallicity and 
a Salpeter IMF, the complete tables for all times, metallicities and 
IMFs will be provided in electronic format. 

The evolution of the accumulated ejecta abundances shows in 
Fig.|4]a sharp increase for Z > 0.008 due to the mass loss of massive 
stars followed by a plateau after the peak of mass loss associated 
with the WR stars decline. For the models with the lowest metal- 
licities, C and O abundances do not change much in relation with 
the initial values, as explained before, decreasing slightly, simulta- 
neously to the increase of N. However, for the other metallicities 
the increase in C and O abundances is quite large. On the other 
hand. He and N abundances, even at the lowest metallicity cases, 
show the sharp increase at 3 Myr followed by a plateau that shows 
an enrichment with respect to the initial value larger than 0.5 dex 
in most cases. This increase in the He and N abundance may have 
important consequences for measurements of the chemical com- 
position of galaxies based on abundances estimated using emission 
lines from Hii regions particularly for the lowest metallicity regions 
i.e. those with 12-l-log(0/H) < 7.3. 

These huge variations in the abundances of the ejecta of a stel- 
lar cluster may also be important for the hydrodynamical evolution 
of the ISM. The high metallicity might lead to extremely short cool- 
ing times in the ejecta with important consequences for the subse- 
quent feedback. 



3 ADDING THE SUPERNOVA EJECTA TO THE 
STELLAR WINDS 

In this section we include the elements produced by supernova in 
our calculations. Stars more massive than in > 12 Mq end their 
evolution as core collapse supemovae. New elements are created 
and ejected in these events. Since we assume an IMF were the most 
massive star m„,„y = 100 Mq has a mean lifetime of t = 3.7 Myr, the 
ejections are zero before this time. Only after 3.7 Myr supernovae 
begin to contribute to the cluster ejected mass. 

In computing the total ejected mass and its composition, an 
important effect to take into account is that due to the stellar wind 
the mass of a star at the pre-supernova stage is smaller than its main 
sequence mass. Thus, the supernova yields for a given star are not 
those corresponding its main sequence mass since the star which 
explodes is less massive. To estimate the supernova yields we took 
as the supernova progenitor mass, the mass of each star at the end 
of its wind phase. In practice the models behave in two different 
ways depending on the wind mass loss rate: 

(i) Small mass loss rate. Stars with initial masses around 15 Mq 
lose only part of their H, therefore X/j > at all times. For example 
stars with Z = 0.02 and initial masses 12, 15 and 20 Mg end the 
winds phase with 11.46, 14.21 and 18.06 Mg, respectively. In this 
case, and even more so for smaller abundances, it is reasonable to 
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Figure 7. Evolution of the instantaneous values of He, C , N, and O abundances in mass in the ejecta of a stellar cluster during the supernova phase. The Une 
codification is as in figure 4. The dotted grey line marks the abundance Z=0.02 in each panel. 



assign them the WW95 models in the same way as PCB98 did. In 
those cases we simply use the final mass of each star to select the 
most appropriate model among the models given by WW95. 

(ii) High mass loss rate. Stars with Z> 0.008 and M> 30 Mohave 
high mass loss rates and arrive to the end of the wind phase without 
H envelope inducing important changes to the supernova explo- 
sio n mechanism and ejecta. For this case we adopted the m odels 
bv lWooslevelail h993h : IWooslev^ Langer & Weaveil ( [1991 here- 
inafter WLW93 and WLW95, respectively). WLW95 calculated the 
evolution of stars between 4 and 20 Mq without H envelope. Their 
results are given for a range of masses of Helium core mnc, defined 
as the mass at which the abundance of H falls to zero. After H ex- 
haustion, these stars continue losing mass. A star of mnc = 20 Mq 
loses 16.44 Mq, ending with a mass of 3.55 Mq. Then, it explodes 
ejecting other 2.00 Mq and keeping 1.55 Mg in the remnant. In 
similar way a star of mnc = 7.00 Mq loses 3.80 Mq ending with 
3.20 Mq before exploding, ejecting 1.70 Mq and producing a rem- 
nant of 1.50 Mq. 

Using the Padova tracks we obtain for the seven original values 
of stellar masses, the mass of each star when H is depleted from the 
envelope. These values are shown in Table |4] for Z> 0.008 given 
that for metallicities lower than this the condition Xh = is never 
reached. The relation between the initial mass in the main sequence 



m, and the mass of the He core is shown in Fig.[5]where we see 
that the behaviour is quite smooth with the resulting ranging 
from 10 to 45 Mg. The explosive yields for masses between 4 
and 20 Mq are taken from WLW95. For masses higher than 20 Mq, 
we take the results from WLW93 where the evolution of WR stars 
is computed starting at Xh=0.50. For models with initial mass of 
60 and 85 Mq, (which begin the RSG or Luminous Blue Variable 
(LBV) phase with 55.5 and 76.9 Mq, respectively), the resulting 
He core masses are mHe= 26.3 Mq and 45.3 Mq, respectively, their 
masses at the end of the wind phase being only 4.25 and 8.30 Mq. 
These models are very similar to our most massive stars of 60 and 
100 Mq, for which m(XH=0.50) are 44.4 and 80 Mq, with mHe = 
27.8 and 44.9 Mq and with final masses of 5.93 and 7.16Mq. 

In Fig. |6] we show the comparison between the evolution 
of the most massive stars used here with the final results of 
WLW93/WLW95. Note that the final masses in the Padova models 
for Z=0.02 are slighfly larger than the ones from WLW95/WNL93 
with the exception of the most massive model. For this case the 
Padova model reaches 7.17 Mq at the end of the evolution while 
the corresponding model from WLW93 has a slightly higher value 
(8.30 Mq). 

Thus, to compute the explosive yields for massive stars with high 
mass loss rate, we proceeded as follows: for each time step of our 
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Figure 8. The evolution of the abundances in mass as log X of the accumulated ejected mass by a stellar cluster during the supernova phase for a) He, b) C, c) 
N and d) O for Z=0.02 and a Salpeter IMF. The dotted grey line marks the reference abundance in each panel. 



Table 5. Accumulated masses ejected by stars dying as supernova explosions for a stellar cluster of 10 
IMF. The complete table with all times, metallicities and IMFs will be provide in electronic format. 



Mo during the first 20 Myr for Z=0.02 and a Salpeter 



IMF 


Z 


time 


m» 


mHc 


mrem 


mej 


H 


He 


C 


N 









(Myr) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


IMF 


0.02 


0.37700E+01 


86.5627 


34.7581 


1.5500 


0.2452E+01 


O.OOOOE+00 


0.1971E+00 


0.4966E+00 


0.3056E-05 


0.1295E+01 


IMF 


0.02 


0.37800E+01 


85.5161 


34.3357 


1.5500 


0.2698E+01 


O.OOOOE+00 


0.2169E+00 


0.5465E+00 


0.3403E-05 


0.1425E+01 


IMF 


0.02 


0.37900E+01 


84.4982 


33.9250 


1.5500 


0.2946E+01 


O.OOOOE+00 


0.2368E+00 


0.5969E+00 


0.3762E-05 


0.1556E+01 


IMF 


0.02 


0.38000E+01 


83.5080 


33.5255 


1.5500 


0.3196E+01 


O.OOOOE+00 


0.2570E+00 


0.6476E+00 


0.4130E-05 


0.1689E+01 


IMF 


0.02 


0.38100E+01 


82.5443 


33.1366 


1.5500 


0.3448E+01 


O.OOOOE+00 


0.2772E+00 


0.6987E+00 


0.4510E-05 


0.1822E+01 


IMF 


0.02 


0.38200E+01 


81.6059 


32.7580 


1.5500 


0.3702E+01 


O.OOOOE+00 


0.2977E+00 


0.7502E+00 


0.4900E-05 


0.1956E+01 


IMF 


0.02 


0.38300E+01 


80.6921 


32.3893 


1.5500 


0.3957E+01 


O.OOOOE+00 


0.3182E+00 


0.8021E+00 


0.5301E-05 


0.2091E+01 


IMF 


0.02 


0.38400E+01 


79.8017 


32.0300 


1.5500 


0.4214E+01 


O.OOOOE+00 


0.3389E+00 


0.8543E+00 


0.5713E-05 


0.2227E+01 



Table|3] we calculate the stellar mass that corresponds to the stellar 
mean-lifetime T(m) = t. Interpolating in table|4]we obtain mnc, i-e. 
the mass of the star at Xh=0. We use this value of mnc to interpo- 
late with the adequate value of mass in the explosive yields from 
WLW95/WLW93 and thus to calculate the stellar yields that corre- 
spond to this star. The yields are then multiplied by the number of 
stars given by 0(m) for the initial (zero age main sequence) mass 
of the star. 

Because the production of elements given by WLW93 and 



WLW95 is calculated only for stars with Z=0.02 we have assumed 
that the relative yields for the other two abundances (Z=0.008 and 
0.05) of the high mass loss rate case are similar. Therefore we have 
calculated from the total ejected masses (the new elements and old 
ones) given by WLW93 and WLW95, the stellar yields, as new el- 
ements ejected by supernova as: 



Pi = meji - (m. - m„,„) * X;,o (3) 
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Figure 9. Evolution of the total abundances of the ejecta of a young stellar cluster with a Salpeter liVIF due to stellar winds and supernova. Stellar wind 
contribution: Blue long dash lines, supernova contribution: Short dashed green lines. Total: Red solid line. The dotted grey line marks the solar abundance in 
each panel. The initial abundances are from top to bottom Z=0,05, 0.02, 0.008, 0.004,0.0004,0.0001 . The elemental abundances are from left to right for He/H, 
C/H, N/H and O/H. 
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where Xi, are the elemental abundances for each element / cor- 
responding to Z=0.02, and m,e,n is the mass of the remnant. 
Then, we recalculated the ejected masses for each Z as: 

meji = Pi + (m. - m„,„) * X, o (4) 

for abundances Z=0.008 and 0.05 using the same stellar yield p, 
and substituting X,_o with the coiTesponding values. 

The instantaneous mass fraction for He, C, N and O for the su- 
pernova ejecta are presented in Fig.|7] Since the ejected mass has no 
H for supemovalb/c, we show all abundances as total mass fraction, 
as before. The dotted line indicates the reference value for Z=0.02. 
Each color indicates a different metallicity with the same coding 
than in previous figures 3 and 4. For the three lowest Z only WW95 
has been used, and so the evolution follows a continuous line. When 
Z > 0.008 we see the change from the ejecta of supernovalb (orig- 
inally WR stars), following WLW93/WLW95 stellar yields, to the 
supemovall from WW95 ejecta as a sharp vertical line. The first 
ones produce less He and N and more C and O than the second 
ones. This behavior is expected given that these stars ejected large 
quantities of He and N before to explode as supernova, during the 
stellar winds phase. 

Table |5] gives the accumulated ejected masses by the super- 
novalb/c and supemovall explosions in each time step. For each 
IMF and metallicity given in columns 1 and 2 we show the time in 
column 3, the mass of the star which dies in this time in column 4, 
the mass of the He core in column 5, the remnant mass in column 
6, the total ejected mass in column 7, and the ejected masses of H 
to O in columns 8 to 12. 

The evolution of the accumulated or total abundances due to 
supernova ejections is illustrated in Fig. |8] A large increase in the 
O and C abundances is clearly seen at 4-6 Myr corresponding to 
the start of the supernova activity. After that there is a steady de- 
crease reaching values around 2 times the initial values after almost 
20 Myr. The He and N abundances, however, do not have an impor- 
tant contribution of supernovalb and they maintain a constant high 
level after about 5 Myr due to the contribution of supemovaH with 
lower mass progenitors. 



4 THE TOTAL ABUNDANCES OF THE CLUSTER 
EJECTA 

In this section we show the evolution of the stellar cluster abun- 
dances obtained by adding both contributions, winds and supernova 
explosions ejections. 

In Fig.[9]we show the time evolution of the abundances. There 
we plot the contribution of winds by long dashed (blue) lines, 
the contribution coming from supernova with short-dashed (green) 
lines and the total abundances with solid (red) lines. For He both 
contributions are more or less similar at the end of 20 Myr when 
Z> 0.004, so the abundance is a factor of two if wind ejections 
are considered compared to the usual calculations performed with 
supernova productions only. For the two lowest metallicities the 
ejected masses are a factor 2 or 3 smaller 

We see that stellar winds produce high abundances of C and 
O only for Z> 0.004. The level of I2+log(0/H) doesn't increase 
for Z< 0.004, while for higher Z it reaches almost two orders of 
magnitude larger. N, however, shows higher abundances than ex- 
pected for all metallicities, even for the two lowest ones. Although 
these results are not unexpected, since they are due to the mass 
loss rate law that depends on Z, these abundances, not calculated 




Figure 10. The dependence of the final ejected masses with the metallicity 
for different IMFs as labelled 

before, may have important consequence over the interpretation of 
observations of Hii regions. 

When we analyze the same plots for supemovae, we see that 
now C and O show very high abundances, compared with the ref- 
erence values, mainly for Z< 0.004, while He and N are roughly in 
the expected level for its metallicity. C and O are primary elements, 
produced directly from He created in the star. Therefore, their pro- 
duction is expected to be independent of the initial composition, 
whereas N is a secondary element produced in the CNO cycle at 
the expense of the initial C and O, so it is reasonable that it scales 
with the initial abundances. 

A summary of these considerations is in Tables|6]and |7]where 
the accumulated masses ejected after 20 Myr by winds and by su- 
pemovae are given. In each one we give for each IMF and each 
metallicity the total mass ejected in Mq and then the contributions 
of each element, H, He, C, N and O proceeding from stellar winds 
and from supemovae, respectively. 

For C the contribution of winds is essential, the supemovae 
contribution being a factor of 10 smaller except for the lowest 
metallicities (Z=0.0001, 0.0004 and 0.004) for which the ejected 
masses are insignificant. Also for N only a small contribution is 
due to supemovae, at the end of the evolution, compared with the 
wind ejecta. For O the supemovae contribution is as important as 
the one from winds, particularly after 10 Myr of evolution. There- 
fore, the contributions of stellar winds to all abundances seems to 
be essential and must be taken into account in the evolution of stel- 
lar clusters and also in the galaxy evolution models in order to in- 
terpret adequately the data. 

We show the final results for different IMFs in Fig. [TO] There 
the fraction of mass ejected in the wind phase in panel a) and the 
final abundances in panels b) to d) are represented as a function 
of Z for the 6 IMFs of this work. Results for MIL and FER show 
the smallest abundances of He, C and N, and the highest H values. 
The other IMFs show results very similar, lower than MIL and FER 
in H, higher for the other elements. In any case differences among 
IMFs results are relatively small. 
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Table 6. Total masses ejected by stellar winds for a cluster at the end of the Table 7. Total masses ejected by supernova for a cluster at the end of the 

wind phase for all abundances and IMFs. supernova phase for each computed metallicity. 



IMF 
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mej 
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He 
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O 


IMF 


Z 


Mej 


H 


He 
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N 


O 






(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 






(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


(Mo) 


SAL 


0.0001 


12787 


9363 


3423 











SAL 


0.0001 


46230 


24090 


17510 


521 


2 


3089 


SAL 


0.0004 


18080 


10881 


7453 


1 


4 


2 


SAL 


0.0004 


46520 


24230 


17670 


529 


2 


3058 


SAL 


0.0040 


41407 


24726 


18598 


17 


71 


47 


SAL 


0.0040 


47510 


24800 


17980 


549 


35 


3018 


SAL 


0.0080 


52781 


27930 


22485 


6933 


160 


2157 


SAL 


0.0080 


49020 


25610 


18690 


563 


72 


2934 


SAL 


0.0200 


62677 


24161 


34006 


9683 


269 


1392 


SAL 


0.0200 


50460 


26510 


19530 


581 


189 


2246 


SAL 


0.0500 


70809 


22410 


40487 


10853 


759 


1446 


SAL 


0.0500 


46240 


25050 


17870 


468 


173 


1522 


MIL 


0.0001 


2862 


2106 


753 











MIL 


0.0001 


36540 


19470 


13850 


410 


1 


2093 


MIL 


0.0004 


4362 


2751 


1646 





1 


1 


MIL 


0.0004 


36860 


19650 


13990 


416 


1 


2076 


MIL 


0.0040 


13721 


8515 


5656 


7 


21 


18 


MIL 


0.0040 


38750 


20630 


14650 


443 


29 


2164 


MIL 


0.0080 


17733 


10159 


6807 


1920 


47 


636 


MIL 


0.0080 


40880 


21760 


15570 


463 


62 


2138 


MIL 


0.0200 


21829 


9598 


10953 


3253 


85 


461 


MIL 


0.0200 


45560 


24180 


17660 


492 


171 


1821 


MIL 


0.0500 


26913 


9890 


13988 


4567 


218 


599 


MIL 


0.0500 


41880 


22880 


16200 


390 


157 


1236 


FER 


0.0001 


2190 


1607 


582 











FER 


0.0001 


16310 


8635 


6179 


183 


1 


978 


FER 


0.0004 


3200 


1968 


1268 





1 





FER 


0.0004 


16440 


8706 


6238 


186 


1 


969 


FER 


0.0040 


8553 


5211 


3679 


4 


14 


10 


FER 


0.0040 


17150 


9078 


6484 


196 


13 


994 


FER 


0.0080 


10966 


6058 


4422 


1316 


31 


422 


FER 


0.0080 


17990 


9531 


6855 


205 


27 


977 


FER 


0.0200 


13278 


5498 


6918 


2021 


54 


288 


FER 


0.0200 


19690 


10420 


7626 


216 


74 


808 


FER 


0.0500 


15719 


5413 


8538 


2556 


146 


337 


FER 


0.0500 


18080 


9860 


6991 


172 


68 


548 


MAR 


0.0001 


32920 


20437 


7472 





1 


1 


MAR 


0.0001 


100900 


52600 


38240 


1138 


4 


6752 


MAR 


0.0004 


46307 


23750 


16270 


2 


8 


4 


MAR 


0.0004 


101600 


52900 


38590 


1156 


4 


6687 


MAR 


0.0040 


98511 


53958 


40587 


37 


155 


102 


MAR 


0.0040 


103700 


54120 


39270 


1198 


77 


6596 


MAR 


0.0080 


128362 


60948 


49078 


15132 


349 


4709 


MAR 


0.0080 


107000 


55890 


40800 


1230 


158 


6410 


MAR 


0.0200 


151042 


52724 


74216 


21128 


587 


3038 


MAR 


0.0200 


110100 


57830 


42610 


1269 


412 


4903 


MAR 


0.0500 


169237 


48903 


88359 


23681 


1656 


3156 


MAR 


0.0500 


100900 


54660 


38990 


1022 


377 


3323 


CHA 


0.0001 


19514 


14286 


5227 





1 


1 


CHA 


0.0001 


66310 


34510 


25120 


748 


3 


4476 


CHA 


0.0004 


27528 


16540 


11381 


1 


5 


3 


CHA 


0.0004 


66720 


34700 


25340 


759 


3 


4430 


CHA 


0.0040 


62281 


37132 


28066 


25 


107 


70 


CHA 


0.0040 


68000 


35440 


25750 


786 


51 


4357 


CHA 


0.0080 


79380 


41850 


33976 


10484 


242 


3254 


CHA 


0.0080 


70060 


36550 


26710 


806 


103 


4231 


CHA 


0.0200 


94122 


36074 


51222 


14535 


406 


2093 


CHA 


0.0200 


71700 


37640 


27750 


830 


269 


3215 


CHA 


0.0500 


105977 


33317 


60835 


16132 


1148 


2155 


CHA 


0.0500 


65690 


35570 


25380 


670 


246 


2179 


KRO 


0.0001 


18392 


13464 


4926 








1 


KRO 


0.0001 


62490 


32520 


23670 


705 


2 


4218 


KRO 


0.0004 


25946 


15589 


10727 


1 


5 


3 


KRO 


0.0004 


62890 


32700 


23880 


716 


3 


4176 


KRO 


0.0040 


58700 


34998 


26452 


24 


101 


66 


KRO 


0.0040 


64090 


33400 


24270 


741 


48 


4107 


KRO 


0.0080 


74816 


39444 


32023 


9881 


228 


3067 


KRO 


0.0080 


66040 


34450 


25180 


760 


97 


3988 


KRO 


0.0200 


88711 


34000 


48277 


13699 


383 


1973 


KRO 


0.0200 


67580 


35470 


26150 


782 


253 


3031 


KRO 


0.0500 


99885 


31402 


57338 


15205 


1082 


2031 


KRO 


0.0500 


61920 


33520 


23920 


631 


232 


2054 


STB 


0.0001 


25645 


18777 


6865 





1 


1 


STB 


0.0001 


92710 


48320 


35120 


1046 


4 


6195 


STB 


0.0004 


36260 


21821 


14947 


1 


7 


4 


STB 


0.0004 


93300 


48600 


35430 


1061 


4 


6133 


STB 


0.0040 


83041 


49587 


37297 


34 


142 


94 


STB 


0.0040 


95280 


49730 


36070 


1100 


71 


6052 


STB 


0.0080 


105851 


56013 


45094 


13903 


321 


4326 


STB 


0.0080 


98310 


51360 


37480 


1130 


145 


5883 


STB 


0.0200 


125698 


48455 


68198 


19418 


540 


2792 


STB 


0.0200 


101200 


53160 


39170 


1165 


379 


4504 


STB 


0.0500 


142007 


44943 


81196 


21766 


1521 


2900 


STB 


0.0500 


92740 


50240 


35840 


939 


347 


3053 



5 CONCLUSIONS 

We have computed the evolution of the total mass ejected and the 
elemental abundances of He, C, N and O for young stellar clusters 
of 10* Mo for 6 IMFs and for 6 different initial metallicities. 

Both stellar wind and supernova contribution to the ejecta are 
included in the computations. The sup ernova contribution include 
the classical gravitational collapse from lWooslev & Weaveij(ll995h 
for low mass progenitors plus the yields calculated by Wooslev et all 
l ll993h and Wooslev, Langer & Weaver ( 1995) for stars more mas- 
sive than 30 Mo and with Z> 0.008 that reach the supernova stage 
after completely losing their H envelope. 

The abundances of the ejecta obtained by adding both contri- 
butions show important differences in comparison with the standard 
method which uses only supemovae yields without taking into ac- 
count the stellar wind yields and/or the difference in the supemovae 



yields resulting from the huge mass loss affecting the pre-supemova 
evolution. 

In particular: 

• The composition of the ejected matter is determined mostly by 
supernova at low metallicities and by stellar winds at around Solar 
metallicities. 

• The total mass ejected by stellar winds ranges from about 1% 
of the initial cluster mass for the lowest metallicity model to about 
6% for the ~ Solar abundance ones of the total mass of the cluster 
for a Salpeter IMF. 

• The total mass ejected by supernova is ~5% of the total mass 
of the cluster for all initial metallicities. 

• At high metallicities the proportion of the mass ejected by the 
winds phase is around 40-60% of the total ejecta. 

• There is a large increase in the abundance of He, C, O and N 
after 2.5 Myr with O and C abundances being the most extreme. 
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The O abundance jumps almost two orders of magnitude between 
2.5 and 4 Myr in our lowest metallicity model and about 3 times for 
the solar abundance model. Between 2 and 3 Myr, the C abundance 
increases between 10 and 30 times its initial value depending on 
the initial abundance. 

• He and N show more moderate jumps than C and O in their 
abundance between 2.5 and 4 Myr. He abundance increases almost 
3 times for the solar value models and about 2 times for Z= 0.0004. 
On the other hand N shows jumps of about 5 times for all abun- 
dances. For cluster ages f < 10 Myr, He and N enrichment is mainly 
due to the stellar winds. 

These huge variations in the abundances of the ejecta of a stel- 
lar cluster can have a profound effect in the hydrodynamical evolu- 
tion of the ISM. The high metallicity of the young cluster ejecta will 
lead to extremely short cooling times with important consequences 
for the subsequent feedback. 

On the other hand, one should keep in mind the main short- 
comings and uncertainties of our models: 

Because present stellar evolutionary models with mass loss do 
not include the supernova phase and supernova models do not cover 
a range in initial abundances, our adopted values for the supernova 
yields of the most massive stars are only approximate. 

Binary evolution is not considered. It is known that in young 
massive clusters perhaps all massive stars are in binary or multiple 
systems, but it is not clear how the presence of a companion would 
affect the properties of wind of a massive star 

The effect of stellar rotation is not included in the stellar evo- 
lution models we have used. Again, as in the case of binaries, is not 
clear at this stage how rotation would affect the wind of a massive 
star. 

In spite of these warnings, our models should be useful for the 
interpretation of the evolution of the ISM in star forming galaxies. 
The resulting tables are available in electronic format. 
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